Nanoradar based on nonlinear dimer nanoantenna 
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We introduce the concept of a nanoradar based on the operation of a nonlinear plasmonic nanoantenna. The 
nanoradar action originates from modulational instability occurred in a dimer nanoantenna consisting of two 
subwavelcngth nonlinear nanoparticlcs. Modulation instability causes a dynamical energy exchange between the 
nanoantenna eigenmodes resulting in periodic scanning of the nanoantenna scattering pattern. Such nanoradar 
demonstrates a wide scanning sector, low operation threshold, and ultrafast time response being potentially 
useful for many applications in nanophotonics circuitry. 
© 2012 Optical Society of America 
OCIS codes: 250.5403, 230.6120, 190.3100, 320.2250. 
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The study of optical nanoantennas is a rapidly devel- 
oping area of research [1]. Nanoantennas keep promises 
for a variety of applications in nanophotonics with ad- 
vantages of enhanced light-matter interaction. For many 
applications, an active control over the nanoantenna ra- 
diation pattern is required. However, the nanoantenna 
operation is usually based on one lasing mode that has 
better directivity and/or higher Purcell factor [2, 3]. 
Nevertheless, for nanoantennas with broken material 
and/or geometrical symmetry, it was suggested to em- 
ploy spectral tunability and variable directionality based 
on switching between different modes [4,5] which also can 
be realized through nonlinear loads [6-9]. 

In this Letter we suggest a novel way for the dy- 
namic control of the nanoantenna directivity by exploit- 
ing modulational instability (MI) of a nonlinear dimer 
nanoantenna. We demonstrate that the development of 
MI in a nanodimer composed of two nonlinear silver 
nanoparticlcs can lead to a periodic exchange of the 
power between the eigenmodes and result in a periodic 
variation of the nanoantenna scattering pattern simi- 
lar to classical phased-array radar systems [10], as il- 
lustrated schematically in Fig. 1. 

We consider a nanoantenna created by a pair of 
nanoparticles placed close to each other and embedded 
into a Si02 host medium with a permittivity Eh = 2.15. 
We assume that the nanoparticles are equivalent, and 
they have the radius a = 10 nm with the center-to- 
center spacing of d = 30 nm (see the insert in Fig. 1). 
Since the condition a/d < 1/3 is satisfied, we can 
employ the dipole approximation [11]. Assuming the 
nanoparticlcs made of silver with a nonlinear Kcrr-likc 
response, we take the dielectric constant in the form 
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\ g + | 2 , where the linear part is given 

by the Drude formula, e\ g = £oo — Wp/[w(w — iv)], with 
4.96, hw p = 9.54 eV, and Uv = 0.055 cV [12], and 
are the local holds inside the first and second par- 
ticles. For spherical silver nanoparticlcs with the radius 



10 nm, we use the nonlinear coefficient with the value 
- 3 x 10 -9 esu from Ref. [13], that is much larger 
than the cubic nonlinearity of SiC-2 (~ 10~ 15 esu [14]), 
the latter is being neglected. Our aim is to study the 
temporal evolution of the scattering pattern of this dimer 
nanoantenna in the nonlinear regime when the response 
of this structure undergoes modulational instability. 

We study a nanoparticlc dimer driven by a plane wave 
(see the insert in Fig. 1) with the frequency close to the 
frequency of the surface plasmon resonance of an indi- 
vidual particle, and analyze the dynamical response of 
the particles through the evolution of their polarizations, 
Pi, 2- To derive the corresponding nonlinear equations for 
the polarizations, we employ the model recently applied 
for the study of strongly localized modes in a chain of 




Fig. 1. (Color online) Schematic of the scattering pat- 
tern of a dimer nanoantenna. Insert shows the structure, 
parameters, and direction of an incident plane wave. 
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nanoparticles [15]. In the case when the external elec- 
tric field E^ ex > is polarized along the dimer axis, this 
model yields the following system of coupled equations 
for the slowly varying amplitudes of the particle dipolc 
moments, 



+ (n - * 7 + \pI I 2 ) Pi + g"p} = e, 
-i^ + {n-ii+\Pi\ 2 )Pi + G»Pl=E, 



(1) 



where the effective coefficient G" = 3£/,./(£ 00 + 
2sh) (a/d) 3 (ikod + 1) exp(— ikod) describes the strength 
of the particle interaction, and we use the dimensionless 
functions 
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and E = -3e fc [x (3) ] 1/2 -B (ex V[ 2 ( £ c» + ^h)?' 2 for the 
slowly varying amplitudes of the particle dipolc mo- 
ments and external electric field, respectively. The in- 
dex Y stands for the longitudinal components with re- 
spect to the dimer axis, the parameter 7 = v/{2ujq) + 
(fcoa) 3 e hi '(e '00 + 2e/j) describes thermal and radiation 
losses of particles, uiq = u> p /(soo + 2sh) 1 ^ 2 is the fre- 
quency of the surface plasmon resonance of an individual 
particle, fco = wo/cy/Eh, £1 = (u> — cjo)/wo, and r = coat. 
For the chosen configuration hu>Q = 3.14 cV. Nonlin- 
ear terms appeared after expressing E-[ 2 via p\ 2 [15]. 
Equation (1) describes the temporal nonlinear dynamics 
of a metallic nanodimer driven by a plane wave with the 
frequency uj ~ loq. 

The stationary states of this system are described 
by homogeneous solutions when the particle dipole mo- 
ments coincide, Pf 2 = Pq ; they can be found as solutions 
of the nonlinear equation, 

H 7 + tt + G" + |P" 



(2) 



where Eq is the stationary amplitude of the external 
electric field. This equation has multiple solutions for 
£1 < —ReG" — y/3 (7 — ImG" ), so that the system of two 
nonlinear nanoparticles is expected to have a bistable 
regime. 

Next, we analyze linear stability of the homogeneous 
stationary states with respect to weak perturbations 
taken in the form of an asymmetric longitudinal eigen- 
mode. By applying a standard technique [16], we derive 
the expression for the instability growth rate, 



A = -ImG" - 7- 



l^'l 4 - 



n - ReG 
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The homogenous state (2) becomes unstable provided 
A(f2, Eq) > 0, and the condition A = defines the bound- 
ary of the MI region on the plane (D,,Eq), as shown 
in Fig. 2. Remarkably, MI occurs for the whole upper 
branch in the bistable region, and it extends much fur- 
ther (see the inset in Fig. 2). 
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Fig. 2. (Color online) Bifurcation diagram on the param- 
eter plane of Eq and fl = (lu—luq)/ujo, with the regions of 
bistability and modulation instability. Inset: Dependence 
of the polarization \Pq\ 2 on Eq at O = —0.1 (along the 
dashed vertical line on the main plot). Red dots mark 
the region of modulation instability. 
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Fig. 3. (Color online) (a) Temporal evolution of the po- 
larizations \Pf\ (red) and \Pj\ (blue), when the exter- 
nal field E (black) approaches and crosses the instabil- 
ity threshold. Instability leads to an energy exchange be- 
tween symmetric and asymmetric eigenmodes. (b) Phase 
trajectories associated with the dynamics shown in (a). 



To study the evolution of this nonlinear system after 
the onset of instability, we perform numerical simula- 
tions of Eq. (1) at zero initial conditions and fixed driving 
frequency fl = —0.12 (hui = 2.76 eV), when the exter- 
nal field E grows slowly approaching the threshold value 
E sa t = 0.014 corresponding to the region of MI. The 
results are presented in Figs. 3(a,b), where we observe 
that instability results in the spontaneous excitation of 
an asymmetric eigenmodc of the dimer leading to an en- 
ergy exchange between the symmetric and asymmetric 
eigenmodes. For weak external fields, the particle po- 
larizations coincide being described by the phase-space 
trajectories corresponding to the focus. When E grows 
reaching the region of MI, this state becomes unstable, 
and the phase trajectories get attracted by two different 
limiting circles, as shown in Fig. 3(b). 

MI breaks the phase locking between the particle 
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Fig. 4. (Color online) Temporal evolution of (a) scat- 
tering cross-section and (b) oscillation angle 9 m i corre- 
sponding to the major scattering lobe, (c-e) Snapshots 
of the scattering intensities (in watts) corresponding to 
maximal deviations from the initial state (red and blue) 
and an intermediate position (green). 



dipolc moments and external field. Consequently, self- 
modulation dynamical response of a nanodimer results 
in temporal modulation of the nanoantenna scattering 
cross-section and rotation of the scattering pattern, as 
shown in Figs. 4(a-e). Period of self-oscillations and con- 
trast of scattering are 90 fs and 18, respectively. This 
modulation is much faster than that typically observed 
using plasmonic antennas with nonlinear semiconductor 
loading [8]. 

We notice a relatively wide angular scanning sector 
6 m i ~ 50° of the nanodimcr's radiation pattern. How- 
ever, using a nanoradar in practice may require much 
better directivity than the directivity of a dimer antenna, 
for which one may utilize, e.g., a high-permittivity dielec- 
tric sphere as a director [2]. 

The saturation amplitude of the external field E sat = 
0.014 corresponds to the intensity of 21.6 MW/cm 2 . Such 
a high illuminating power can lead to thermal damage 
unless we operate with short pulses. To estimate the 
maximal pulse duration, we use the value of the abla- 
tion threshold of 3.96 J/cm 2 obtained for silver particles 
in a Si02 host matrix in the picosecond regime of illumi- 
nation [17]. Taking into account the amplification of the 
electric field inside the Ag nanoparticle due to surface 
plasmon resonance and the required intensity, we evalu- 
ate the maximal pulse duration as 0.5 ns, which is much 
longer than the characteristic period of self-oscillations. 
Thus, the predicted scattering pattern rotation can be 
readily observed in experiment. 



In conclusion, we have suggested a novel approach for 
realizing an active control over the optical nanoantenna 
directivity. We have shown that modulational instability 
in a nonlinear plasmonic dimer can lead to a dynamical 
energy exchange between the nanoantenna eigenmodes 
resulting in a periodical variation of the nanoantenna 
scattering pattern. Such actively tunable nanoantennas 
may operate in a wide scanning sector with a low op- 
eration threshold and ultrafast modulation. Further op- 
timization may be possible by considering other designs 
of nanoantennas with higher directionality, such as Yagi- 
Uda [1] or metal-dielectric [2] configurations. Another 
possibility for tunability of the nanoradar is a variation 
in the plane wave incidence, intensity, and frequency. In 
this case, different scenarios of the development of mod- 
ulational instability include transitions to other steady 
states, bifurcations, period doubling, and even a transi- 
tion to chaos. We have studied some of such regimes and 
results will be published elsewhere. 
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